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Abstract 


The generalized mathematical formulation of sloshing dynamics for partially 
filled liquid of cryogenic superfluid helium II in dewar containers driven by the 
gravity gradient and jitter accelerations associated with slew motion for the 
purpose to perforin scientific observation during the normal spacecraft operation 
are investigated. An example is given with the Advanced X-Ray Astrophysics 
Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the 
sloshing dynamics. The jitter accelerations include slew motion, spinning 
motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising 
from machinery vibrations, thruster firing, pointing control of spacecraft, crew 
motion, etc. Explicit mathematical expressions to cover these forces acting on 
the spacecraft fluid systems are derived. The numerical computation of sloshing 
dynamics is based on the non- inertia frame spacecraft bound coordinate, and solve 
time -dependent, three-dimensional formulations of partial differential equations 
subject to initial and boundary conditions. The explicit mathematical 
expressions of boundary conditions to cover capillary force effect on the liquid- 
vapor interface in microgravity environments are also derived. The formulations 
of fluid moment and angular moment fluctuations in fluid profiles induced by the 
sloshing dynamics, together with fluid stress and moment fluctuations exerted on 
the spacecraft dewar containers have also been derived. Examples are also given 
for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with 
slew motion. 
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Statement of Work for the Numerical Simulation of 
The AXAF-S Superfluid Helium Cryogenic Dewar Slosh Dynamics 

In order to support the Advanced X-ray Astrophysics Facility- Spectroscopy 
(AXAF-S) Cryogenic Subsystem (CSS) Preliminary Design Review (PDR) , the following 
tasks have been carried out and accomplished in this research project: 

1. Adopt the numerical model of the Gravity Probe -B cryogenic, 
superfluid helium dewar to the AXAF-S configuration by removing the 
center core tube, adjusting fluid dynamic coefficients for AXAF-S 
nominal fluid temperature of 1.3 K, and adjusting the rigid body 
spacecraft mass properties and orbital parameters to those associated 
with AXAF-S. In compliance with AXAF-S numerical model, mathematical 
formulations of the dynamic model of gravity gradient and external 
jitter accelerations associated with slew motions have been 
developed. Detailed descriptions are illustrated in pages 5 to 27 
of the present report. 

2. Assume a worst case ratio of fluid- to -dewar volume with respect to 
slosh dynamics in the range of 70% and establish the equilibrium 
conditions for fluid state with zero vehicle internal angular 
velocity. Figure 10 in this report shows the equilibrium conditions 
illustrated . 

3. In order to perform the numerical simulations to determine fluid 
profiles, the levels of force and torque applied to the AXAF-S 
spacecraft by the superfluid in response to the gravity gradient and 
jitter accelerations associated with slew motion. Mathematical 
derivation of fluid profiles including angular momentum and moment 
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fluctuations driven by both gravity gradient and j itter accelerations 
associated with slew motions are illustrated in pages 27 to 30 while 
the mathematical formulation of force and torque applied to the 
spacecraft system are described in pages 30 to 37 of the present 
report. Examples of numerical simulations applicable to AXAF-S 
spacecraft by superfluid helium in response to vehicle accelerations 
have been carried out which include the following items: 

(a) Time -dependent, three-dimensional bubble configuration 
oscillations, bubble mass center fluctuations, angular momentum and 
moment fluctuations of 70% liquid-f illed superfluid helium flow 
profiles, and force and torque actuated on the dewar container of 
spacecraft driven by gravity gradient acceleration associated with 
slew motion during a time period of 800 seconds have been carried 
out. The results are fully illustrated in pages 37 to 53 and in 
Figures 10 to 15 and 22 to 27 of this report. 

(b) Time - dependent , three-dimensional liquid-vapor interface 
oscillations of superfluid helium, bubble mass center fluctuations, 
angular momentum and moment disturbances of 70% liquid- filled 
cryogenic helium flow profiles, and force and torque applied on the 
dewar container of AXAF-S spacecraft driven by external acceleration 
jitter associated with slew motion at amplitudes of 10’ 5 ^ (g 0 -9.81 
m/s 2 ) during a time period of 800 seconds have been carried out. 
These results have been illustrated in pages 39 to 53 and in Figures 
16 to 21 and 28 to 31 of this report. 
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Because unexpected complications in numerical instabilities and convergent 
problems occurred during the execution of iterations in numerical computation, 
computation of white noise frequency spectrum was unable to accomplish during the 
period of time in this research. The total expenses of computer service are 
completely covered by the University of Alabama in Huntsville (UAH) without 
charge to the government. The amount of computational expenses supported by UAH 
exceeded the amount of expenses supported by the government in this project. 

Data outputs from these numerical simulations had been provided to the 
Contracting Officer's Technical Representative (COTR) in both graphical form and 
on 3 . 5 inch micro -floppy diskettes in a disk format as required by the contract. 
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I. Introduction 


For the purpose to carry out scientific experiments, some experimental 
spacecraft use cryogenic cooling for observation instrumentation and telescope, 
superconducting sensors for gyro read-out and maintain very low temperature near 
absolute zero for mechanical stability. The approaches to both cooling and 
control involve the use of superfluid liquid helium II. In this study, sloshing 
dynamics associated with spinning and/or slew motions are investigated. To cover 
the spacecraft spinning and/or slew motions, Gravity Probe-B (GP-B) and Advanced 
X-Ray Astrophysics Facility- Spectroscopy (AXAF-S) spacecrafts have been chosen 
as the examples in this study. Both the GP-B and the AXAF-S spacecrafts adopt 
the cooling and boil-off from the cryogenic liquid helium dewar as a cryogen and 
propellant to maintain the cooling of instrumentations, attitude control and 
drag-free operation of the spacecraft. The potential problems for cryogenic 
liquid in dewar container could be due to asymmetry in the static liquid helium 
distribution and to perturbations in the liquid-vapor interface caused by slosh 
wave excitation driven by pointing control, machinery vibration, etc. 

For the cases of both the GP-B and the AXAF-S spacecrafts, cryogenic liquid 
helium II, at a temperature of 1.3 K, is used as the propellant. With its 
superconducting behavior, there is no temperature gradients in the liquid helium. 
In the absence of temperature gradient along the surface which drive Marangoni 
convection, the equilibrium shape of the free surface is governed by a balance 
of capillary, centrifugal and gravitational forces. Determination of liquid- 
vapor interface profiles based on computational experiments can uncover details 
of the flow which can not be easily visualized or measured experimentally in a 
microgravity environment. 

The instability of the liquid-vapor interface surface can be induced by the 
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presence of longitudinal and lateral accelerations. Slosh waves are, thus, 
excited which produces high and low frequency oscillations in the liquid 
propellant. The sources of the residual accelerations range from the effects of 
the Earth's gravity gradient and jitter accelerations which include, atmospheric 
drag on the spacecraft, vibration of compressor, spacecraft attitude motions 
arising from machinery vibrations, thruster firings, spacecraft slew motion, 
pointing control of spacecraft, crew motion, etc. Recent study (Kamotani et 
al., 1981) suggest that the high frequency accelerations may be unimportant in 
comparison to the residual motions caused by low frequency accelerations. 

Time-dependent dynamical behavior of partially-f illed rotating fluids in 
reduced gravity environments was simulated by numerically solving the Navier 
Stokes equations subject to the initial and the boundary conditions (Hung and 
Shyu, 1991 a ,b , c ; 1992 a,b,c; Hung et al., 1991 a,b,c; 1992 a,b,c). At the 
interface between the liquid and the gaseous fluids , both the kinematic surface 
boundary condition, and the interface stress conditions for components tangential 
and normal to the interface, were applied (Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 
Hung et al., 1991 a,b,c; 1992 a,b,c). The initial conditions were adopted from 
the steady- state formulations developed by Hung et al (1989 a,b,c) . Some of the 
steady- state formulations of interface shapes were compared with the available 
experiments carried out by Leslie (1985) in a free-falling aircraft (KC-135). 
The experiments carried out by Mason et al (1978) showed that the classical fluid 
mechanics theory is applicable for cryogenic liquid helium in large containers. 

In the spacecraft orbit around the Earth, the direction of azimuth angle 
of Earth toward the location of the spacecraft geometric center varies from 0° 
along the rolling axis of spacecraft to various directions in which three 
dimensional calculation shall be adopted. 
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As the spacecraft moves along the orbit, any fluid capable of motion 
relative to the spacecraft is subject to the acceleration that arises from the 
gravity gradients of the Earth (Avduyevsky, 1984; Forward, 1982; Misner et al . , 
1973). The interaction between the particle mass of fluid and the spacecraft 
mass due to gravity gradient accelerations (Forward, 1982) are capable for the 
excitation of slosh waves and disturb the fluid system which induces the 
fluctuations of viscous stress and its moment exerted on the containers of the 
spacecraft fluid systems. In the meanwhile, the sources of residual acceleration 
of gravity jitter range from atmospheric drag on the spacecraft, background 
gravity, spacecraft attitude motions arising from machinery vibrations, 
spacecraft slew motion, thruster firings, crew motion, etc., are also capable for 
the excitation of slosh waves on the fluid containers. 

It is critically important to understand the physical and dynamical 
behavior of cryogenic helium in a rotating cylinder to effectively promote space- 
oriented missions. 

At temperatures close to absolute zero, quantum effects begin to be of 
importance in the properties of fluids. At a temperature of 2.17°K, liquid 
helium has a A-point (a second-order phase transition) ; at temperatures below 
this point liquid helium (helium II) has a number of remarkable properties, the 
most important of which is superfluidity. This is the property of being able to 
flow without viscosity in narrow capillaries or gaps. 

The basis of the dynamics of helium II is the following fundamental result 
of microscopic theory. At temperatures other than zero, helium II behaves as if 
it were a mixture of two different liquids. One of these is a superfluid and 
moves with zero viscosity along a solid surface. The other is a normal viscous 
fluid. It is of great importance that no friction occurs between these two parts 
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of the liquid in their relative motion, i.e. , no momentum is transferred from one 
to the other. 

It should be emphasized that regarding the liquid as a mixture of normal 
and superfluid parts is no more than a convenient description of the phenomena 
which occur in a fluid where quantum effects are important. One of these motions 
is normal and has the same properties as the motion of an ordinary viscous fluid, 
but the other is the motion of a superfluid. The two motions occur without any 
transfer of momentum from one to another. We can, in a certain sense, speak of 
the superfluid and normal parts of the fluid, but this does not mean that the 
fluid can actually be separated into two such parts. 

If fluid flow can separate helium II into the regions of the superfluid and 
the normal fluid, two temperatures zones are immediately created. A very low 
temperature zone is located at the zone of very high density concentration of the 
superfluid, while a high temperature (below 2.17° K) zone is located at the zone 
of very high density concentration of the normal fluid at the other end. The 
existence of a sharp temperature gradient at the interface between the superfluid 
and the normal fluid results in the creation of a great difference in chemical 
potential, which, in turn, induces a great reverse pressure gradient, creating 
the environment of isothermal fluid distribution everywhere throughout the 
cylinder and homogenous distribution of superfluid density concentration. This 
illustration of the possible separation of superfluid from normal fluid of helium 
II means that there is in reality, no way for anyone to achieve the separation 
of the superfluid from the normal fluid of helium II. In other words, in 
considering the dynamical behavior of helium II in a large rotating cylinder, a 
mixture of the superfluid and the normal fluid without separation of the two 
fluids is accounted for in the model computation. Density concentration of 
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superfluid is a function of temperature, which is also true for the surface 
tension and viscous coefficient for helium II (Wilks, 1967; Hoare et al. , 1961; 
Hung, 1990; Hung and Lee, 1992). In fact, the experiments carried out by Mason 
et al. (1978) showed that the classical fluid mechanics theory is applicable for 
cryogenic liquid helium in a large container. In this study, the theory of 
viscous Newtonian fluids is employed with modification of transport coefficients 
adjusted by normal and superfluid density concentration which is a function of 
temperature . 

II. Functions of Scientific Observation and Spacecraft Motions 

The AXAF-S spacecraft is a sun synchronous Earth satellites orbiting at 650 
km altitude directly over the poles. The functions of scientific observation for 
the AXAF-S spacecrafts and its motions are illustrated as follows: 

The AXAF-S with its sister spacecraft AXAF-I (I for Imaging) are two 
spacecrafts restructured from the original AXAF design to carry out astrophysical 
observation. Equipped with the (microcalorimeter) X-Ray Spectrometer (XRS), 
AXAF-S provides high-throughput, high-resolution, non-dispersive spectroscopy at 
the higher AXAF x-ray energies - including the astrophysically important iron-K 
spectral region (above 6.4 keV) - and also permits some spatially resolved high- 
resolution spectroscopy. AXAF-S comprises a foil-mirror (or possibly, a 
replication-optic) telescope (4.7-m focal length), with XRS in the focal plane. 
With the baseline optical system, the AXAF-S provides important, unique 
capabilities for high-throughput, high resolution (above 1 keV) spectroscopy of 
extended and point sources and for some spatially resolved high- resolution 
spectroscopy . 

Because of the unique capabilities for high-resolution spectroscopy of 
point and extended sources, AXAF-S is eligible to carry out following scientific 
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programs : 


(1) Distance scale: Measure the cosmological distance scale using the 

Krolik-Raymond-Sarazin technique in clusters of galaxies. The 
technique requires comparison of the equivalent width of an 
absorption line of the spectrum of a background (or central) active 
galactic nuclei (AGN) observed through a cluster of galaxies, with 
emission- line flux of the same line from the cluster itself. 

(2) Active Galactic Nuclei (AGN): Detect x-ray emission from the hot 

intercloud medium around AGN. Study photo- ionized environs of AGN 
through analysis of the redshift and line shape of the resonance 
fluorescence line of iron. 

(3) Absorption features: Analyze absorption edges to determine 

abundance, ionization, and column densities of absorbing clouds 
around active galactic nuclei and in the interstellar medium. 

(4) Cooling flows: Conduct a detailed study of cooling flows in giant 

elliptical galaxies and in clusters of galaxies. Isolate emission 
from cooling regions from the hot plasma. Determine abundances, 
temperatures, and their radial dependence. 

(5) Clusters of galaxies: Measure redshift and velocity dispersion of 

intracluster medium. Determine the dynamical properties of the 
emitting gas, for direct comparison with simulations of cluster 
merging. Measure the evolution of the abundances of low-Z elements. 

(6) Supernova remnants: Trace plasma diagnostics in low-surface- 

brightness areas of supernova remnants. Determine spatial dependence 
of temperature and ionization age for young supernova remnants, using 
spatially resolved high-resolution spectra. 
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(7) X-ray binaries: Determine the sites at which Fe line emission 

originates in individual binary x-ray sources. The width and 
temporal behavior of (hot-plasma) recombination and (cold-plasma) 
fluorescence lines (most strongly Fe-K) probe the spatial structure 
of the environs of the compact x-ray source. 

(8) Stars: Perform time-resolved spectroscopy to examine plasma heating 

and cooling during stellar flares. For eclipsing binaries, determine 
location of active regions of the corona. 

To comprise these functions of scientific observation, AXAF-S equipped with 
foil-mirror telescope and XRS in the focal plane is capable to observe point and 
extended sources of AGN, clusters of galaxies, supernova remnants, x-ray 
binaries, stellar flares, active regions of corona, etc., through spacecraft slew 
motion of pointing control. In other words, spacecraft slew motion without 
spinning with rotating axis is required for AXAF-S to perform its scientific 
mission. 

III. Basic Characteristics of Gravity Gradient 
and Gravity Jitter Accelerations 

Any fluid element inside the on-orbit spacecraft fluid system is subject 
to the acceleration that arises from the gravity gradient of the Earth 
(Avduyevsky , 1984; Forward, 1982; Misner et al., 1973; Hung and Pan, 1993; Hung 
et al., 1993 a,b,c). Once the spacecraft orbit is fixed, the orbit period is 
determined and the basic structure of the gravity gradient acceleration also can 
be calculated. However, gravity gradient acceleration acting on each fluid 
element inside the on-orbit spacecraft fluid system is different dependent upon 
the distance of the location of the fluid element to the geometrical center of 
the spacecraft and its direction toward the location of the center of the Earth. 
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This acceleration can only be calculated based on the non- inertia frame of 
spacecraft bound coordinate. Thus, the coordinate system shall be transformed 
from ordinary inertia frame coordinate to non- inertia coordinate. 

(A) Orbit Motion of Spacecraft 

Let us consider the case of the AXAF-S spacecraft, which is the Earth 
satellites orbiting at 650 km altitude directly over the poles, the orbit period, 
r Q can be computed from following expression: 


t o =2* 


rV 2 

R E gV 2 


(3-1) 


where R E denotes the radius of Earth (- 6373 km); R c , the radius of the circular 
orbit (■* Rg + h — 7023 km); h, orbit altitude (— 650 km); and gQ, Earth gravity 
acceleration (- 9.81 m /s z ) . For the case of the AXAF-S spacecraft, the orbit 
period r 0 - 97.6 min, and orbit rate n - “ 1-07 x 10“ 3 rad / a . 

As the spacecraft is orbiting around the Earth, the azimuth angle of the 
Earth, toward the location of the spacecraft geometric center varies with 
respect to time. At time t - 0, the rolling axis of the spacecraft is aligned 
with the radial direction of the Earth's center to the spacecraft geometric 
center. Assuming the spacecraft rolling axis is linearly turning around 0° to 
360° in the orbit period, r c , of the spacecraft when the spacecraft is orbiting 
around the Earth. Without the spacecraft slew motion, the azimuth angle (v^e 0 ) 
can be defined as 




2tc 


(3-2) 


where r Q is the spacecraft orbit period [defined in Equation (3-1)] ; and t is the 
time measured from the instant when the direction of the spacecraft rolling axis 
is aligned with the radial direction of the spacecraft mass center to the center 
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of the Earth. 


(B) Slew Motion of Spacecraft 

For the purpose to carry out wide-range observations, some scientific 
spacecraft requires slew motion with respect to the mass center of the 
spacecraft. This is particularly true for the case of the AXAF-S spacecraft. 
For the case of the spacecraft slew motion, azimuth angle, shown in Equation (3- 
2) , shall be modified through the coordinate transformation of slew motion when 
the spacecraft is orbiting around the Earth. 

Let us assume that the slew motion starts with the center located at the 
mass center of the spacecraft. Let us choose cartesian coordinate (x" , y" , z") 
with z"-axis along the axis of the dewar container (see Figure 1). At time t - 
0, the radial vector r c from the center of the spacecraft to the center of the 
Earth lies on the x M -z" plane of the cartesian coordinate chosen (see Figure 2). 
The azimuth angle ^ is defined as the angle between the radial vector r c and the 
z"-axis. Rotation matrices for spinning and/or slew motions along the x"-, y"- 
and z"-axes can be expressed as 


1 

0 

0 


COSO) y t 

0 

-siruOyt 


COS6) z t 

sino> z t 

O' 

0 

COSO) x t 

sinto^t 


0 

1 

0 


-sin 

c osu z t 

0 

0 

-sino^fc 

COS(i)„t 
x - 


sina> y t 

0 

COSGJyt 


0 

0 

1 


respectively. Here, o> x , w y and u> z denote angular velocity of slew and/or 
spinning motions along the x"-, y M - and z"-axes, respectively. Radial vector r c 
in cartesian coordinate without slew and spinning motion is (see Figure 2) 

f co = [sim|f £o< 0, -cosi|f Fo ] (3-3) 

With an execution of spinning motion along the z"-axis only, radial vector r c 
becomes (see Figure 2) 
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cos<i) z t sino) r t 0 


_ 



sini|/ £o 

-sino> z t cosu z t 0 


0 

0 0 1 


-cos^ 


= [simjr^cosc^t, -sinijf^sino^t, -cosi|r So ] (3-4) 


With an execution of slew motion along the y"-axis only, radial vector r c becomes 
(see Figure 2) 


■c-y 


cosu> y t 0 -sino> y t 


r 



sinijr*. 

0 10 


0 

sinw y t 0 cosco y t 


-cosily 


= fsin(ijf so + to y t) , 0, -cos (ijr So + a) y t) ] (3-5) 

With an operation of slew motion along the x"-axis only, radial vector r c becomes 
(see Figure 2) 


10 0 
0 cosu^t sinG> x t 
0 -sinoj^t cos(.) x ti 


sini|f eo 

0 


-cosi|» B ! 


= [sinijr^, -cosijr^sinoj^t, -cosilr^cosu^t] (3-6) 

In other words, radial vector r c will be modified from the mathematical 
expression shown in Equation (3-3) to (3-4), (3-5) and (3-6) for the slew and/or 
spinning motions along the z" -, y"- f and x" -axes alone, respectively. In 
particular, for the case of slew motion along the y"-axis, comparison between 
Equations (3-3) and (3-5), it shows that the azimuth angle will be modified as 

+ (3-7) 

For the successive operations of the spacecraft from spinning motion along 
the z"-axis, then slew motion along the y"-axis, and then slew motion along the 
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x"-axis, radial vector r c results 


■ c-z,y,x 


1 

0 

0 


0 

1 


COS(i> y t 

0 

-sin.<o y fc 

COSCD^t 

sina) x t 


0 

1 

0 

-sino> x C 

coso x t 


sinw y t 

0 

COS(O y t 


coso) z t sinw r t 0 
• -sinti ) r t coso) z 0 
0 0 1 


si.in|r So 

0 


-cosy Eo 


(3-8) 


To give a practical example for spacecraft slew motion along the x"-axis 
and/or y"-axis, starting and ending conditions for angular position, angular 
velocity and angular acceleration of slew motion can be given as 0 s (t o ) - 0, 
*.(t f ) - */2, MO - M t f ) - t % (t 0 ) - H t f ) - tf(t 0 ) - T( t f ) - 0 where t 0> t f , 

f4i 

Os, 0 S , and denote starting time, ending time, first time derivative, second 
time derivative and third time derivative of the angle of slew motion, 
respectively, expressions of 0 a become 


0 5 (tr) = a 1 t 1 + a 6 t 6 + a 5 t s + a 4 t 4 


(rad) (3-9) 


6 3 (t) = Q s (t) 


la^t 6 + 6a 6 t 5 + 5a s t 4 + 4a 4 t 3 (rad/ min) (3*10) 


rl 

0 s (t) = a s (t) = 42 a 7 t 5 + 30a 6 t 4 + 20a s t 3 + 12a 4 t 2 (rad/min 2 ) (3-11) 

where o> s and a 5 denote angular velocity and angular acceleration of slew motion, 
respectively. If the slew motion operates at 90° in 10 minutes, the coefficients 
of angular position a 7 , a 6 , a 5 and a 4 , shown in Equations (3-4) to (3-6), can be 
expressed as follows: 

a 7 - -3.14159 x 10‘ 6 , a 6 - 1.09956 x 10“\ a 5 - -1.31947 x 10" 3 

a A - 5.49778 x 10" 3 , and t f - 10 min. 

The expressions of 0 a , u> s , a s are illustrated in Figures 3, 4 and 5, 
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respectively . 

In addition to the modification of the azimuth angle made by the spacecraft 
slew motion through the formulation of coordinate transformation, shown in 
Equation (3-3) to (3-8), accelerations are also induced to activate on the fluid 
mass in the dewar container. Accelerations acting on the fluid particle in the 
dewar induced by the slew motion of the spacecraft with the coordinate fixed at 
the spacecraft center of the mass is as follows (see Figure 1): 

«» 

R p = (3 x(G>xR p ) + axR p + 2G>xv (3-12) 

where ^ denotes the position vector of the fluid particle in the dewar container 
relative to the body frame of the spacecraft ; "cT, angular velocity of the 
spacecraft body frame; cT, angular acceleration of the spacecraft body frame; and 
v> velocity of the fluid particle relative to the spacecraft body frame. 

As we indicated earlier, let us assume that the slew motion starts with the 
center located at the spacecraft mass center, cartesian coordinate (x' r y' ' 
z'') is chosen with origin located at the spacecraft mass center. Let us also 
assume that plane intersects the center of Earth and the spacecraft mass 

^®nter. In other words, azimuth angle of Earth toward the spacecraft mass center 
lies in the x''-z'' plane. Slew motion is along both the x' ' - and y''- 

^ m » 

coordinates. Thus, w s - (w sx , cj sy , 0) and a s - (a sx , c* sy , 0) , due to slew motion 
becomes 


#4 

R x " 



R p,slew ~ 

Ryll 

= 

"^x W jy) ~ a sx^z ~ ^ ^sx^z 


R z n 

slew 

- R z^J+u sy 2 ) + (a^ y - a sy R x ) + 2(u sx v y -o sy v x ) 


( 3 - 13 ) 

(C) Coupling for the Accelerations of Spinning and Slew Motion of Spacecraft 
As we indicated in Section II of this report, with their specified 
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functions of scientific observation, the dewar container of the GP-B is spinning 
with a certain rotating rate without slew motion during the normal operation 
while AXAF-S requires slew motion for pointing control to observe point and 
extended sources of astronomical objects without spinning. For some particular 
reasons required in other spacecraft, it might be faced with the situation that 
both spinning and slew motions are needed simultaneously. To encounter this 
case, the following formulations are made to deal with coupling for the 
accelerations of spinning and slew motion of the spacecraft: 


slew end spinning 


slew and spinning 


u sy {u> sx Ry-u sy R x ) + <x sy R z + 2 (o sy v z 

^ sx ( ^ sx^y*^ sy^x^ ~ ® sx^ z ~ sx V z 

+ »«K a > + ( a sx R y~ (i sy R x > + 2 (' ^ sx V y " sy V X ) 


slew 


-2a> z v y 

■2b) z v x (3-14) 

- spinning and coupling 

where and u> z denote angular velocity and angular acceleration, respectively, 
of spacecraft spinning motion along the z-axis. 

For the case of the GP-B spacecraft, there is no slew motion and the 
spinning is the only acceleration acting on the spacecraft fluid system. 
Acceleration due to spacecraft spinning motion becomes 


< “sy* z - W z*y> <*>z + <•> J^x * 




0 ). 


p, spinning 


Rx 


- R x<* z 2 -R y u z - 2(O z V y 

Ry 

= 

-R y a 2 2 +R x u z + 2u z v x 


spinning 

0 



(3-15) 


Ispinning 

To convert the expression of Equation (3-15) in cartesian coordinate to 
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cylindrical coordinate, by using the relationships of (R,, Ry) - (rcostf, rsin*) 
and (v x , v y ) - (u r cos0 - u e sin0 , u r sin0 + u 9 cosd) t Equation (3-15) becomes 


** 

R x 


-rcos0o>/-rsin0d) z -2 (u r sin0 + u s cosQ) u z 

JD 

p, spinning ~ 

Ry 

= 

-rsin0Q z 2 +rcos0d) z + 2 (u r cos0 - i^sin.0) <j z 


R x \ 

spinning 

0 


and (3-16) 


R — 

* x p, spinning “ 

*r 

R* 

— 

* m m 
J? x cos0 + i? y sin0 

-J? x sin0 + i? y cos0 

. 

-ro) z 2 -2u a G) z 
rw g + 2u r w z 

(3-17) 


R z 

spinning 

R z 

spinning 

0 

spinning 


Accelerations induced by spacecraft spinning motion alone becomes 


’ a r 


Rr 


r<o z 2 +2UeO> z 

a e 

= — 

R« 

= 

-ru z -2u r « z 

a z. 

spinning 

sj 

spinning 

0 


(3-18) 


spinning 


(D) Gravity Gradient Acceleration 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 


d gg =-n 2 [2 (f c ‘d)f c -d] (3-19) 

where a gg denotes gravity gradient acceleration vector; d, the vector (not a unit 
vector) from the fluid element to the spacecraft geometric center; r c , a unit 
vector from the spacecraft geometric center to the center of the Earth; and n, 
the orbit rate (see Figure 2). 

It is assumed that the gravity gradient exerted on the geometrical center 
of the spacecraft orbiting around the Earth on its specified orbit is zero. In 
other words, all the gravity acceleration exerted on the spacecraft is nothing 
but the gravity gradient acceleration which is defined in Equation (3-19). In 
this study, we are interested in investigating how gravity gradient acceleration 
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affects the dynamical behaviors of cryogenic fluid elements of helium. 

For the convenience of mathematical calculation, let us describe all the 
parameters involved in Equation (3-19) in terms of cartesian coordinates. In 
order to match with the computer simulation, mathematical derivation are 
considered in the first quadrant. Figure 2 illustrates the geometrical 
relationship of the parameters shown in Equation (3-19). 

Let us consider the fluid element of interests, m, located at (r, 6, z) in 
cylindrical coordinates and at (x, y, z) in cartesian coordinates. The origin 
of the two coordinate systems is located at the center bottom of the dewar tank. 
The slew and/or spinning motions, mentioned earlier, are executed at the 
spacecraft mass center with cartesian coordinate (x" , y" , z"). The geometry 
center or the spacecraft mass center is located at z - L c . As |d| (not an unit 
vector) is much smaller than the distance between the location of the AXAF-S 
spacecraft geometric center to the center of the Earth, r c (an unit vector) 
through the AXAF-S geometric center and r e (an unit vector) through the fluid 
element, m, is basically the same. Assume that vector r c lies in the x-z plane 
of the cartesian coordinate. 

Radial vector r c with the modification of slew and/or spinning motions 
along the x"-, y"-, z"-axes have been derived in Equations (3-3) to (3-8). Based 
on the relationship between coordinates (x, y, z) and (x" , y" , z") 

i fl 0 0] x /i 

= 010 y" (3-20) 

0 0 lj \z". 

Vector & in (x, y, z) coordinate becomes 

d = [-rcosO, -zsin0, - (z-L c ) ] (3-21) 

Substituting Equations (3-4) and (3-21) in (3-19), non-inertia frame 


X 

y 

z-L 
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expression of gravity gradient acceleration with spinning motion in z-axis 


becomes 


a 

a 

* 


99. r | 

99. * 


spinning la x-*jcix 


-n 2 


3 [ -rsinf ^cos (0 + <*> f t) ♦ ( z~L c ) cos^^] sin^r^coaw^t ♦ rcos8 
-3 [-rsint^cos (8 + w 4 t) ♦ (r-L e ) cost *,] ®ini|r Jto Hintt x t *rsin0 
-3 t-rsinilr^cOB (0 + « x t) ♦ {z-L c ) costj*] cost*,* (z-L c ) 


(3-22) 


Substituting Equations (3-5) and (3-21) in (3-19), non-inertia frame 
expressions of gravity gradient acceleration with slew motion in y-axis becomes 


a 99.x 


3 [-rsini|rcos0 +cosi|r ( z-L c ) ] sirujr +rcos0 

a 99.y 

= —23 2 

rsin0 

^99. 2. 

slew in y-axis 

-3 [ -rsini|fcos0 + cosi|f (z-L c ) ] cost|r + ( z-L c ) 


(3-23) 

where ^ + « y t. 

Substituting Equations (3-6) and (3-21) in (3-19), non-inertia frame 
expressions of gravity gradient acceleration with slew motion in x-axis becomes 


a 

a 


99.x 
99. Y 


1*2 •» In jc-ajcIx 


- n 2 


3{-rcos0sin\^ JO >cos^ ab [rsin<i>^tsin0+cosw x t ] }sinf Jto +rcoB8 

*-3{-rcos0flint*,>cos\jr to ) ] ain« x tsin0+cosM x t (z-L e ) ] Tcost to sinu x £: ursine 
-3{-rcos0sini|r ft> >cost la [rsin(«> Jt tsin0-«'COSw x t (z-Z. c ) ] }co8t fc coso> x t+ (z-L c ) 


(3-24) 


The relationship for the coordinate transformation from cartesian to 
cylindrical coordinates for any vector F (such as velocity or force vectors ) in 
non- inertia frame of spacecraft bound coordinate can be shown as 


\F r ] 


cos0 sin0 0 


pwi 


=S 

-sin0 cos0 0 


F y 

F *\ 


0 0 1 




Thus, the gravity gradient acceleration located at (r, $ , z) can be computed from 
that located at (x, y, z) , shown in Equations (3-22) to (3-24), from the 
following relation: 



_1 


cos0 sin© 0 


3L 


99, r 

a 99.& 


-sin0 cos0 0 


99, x 

a 99.y 


a 99. 2. 


0 0 1 


a 99. 2. 


(3-26) 
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(E) Jitter Accelerations 


In addition to gravity gradient acceleration acting on the fluid element 
of on-orbit spacecraft fluid systems, there is another acceleration of gravity 
jitter also exerted forces on the fluid systems. The sources of residual 
acceleration of gravity jitter range from slew motion of spacecraft, atmospheric 
drag on the spacecraft, background gravity, spacecraft attitude motions arising 
from machinery vibrations , thruster firings, crew motion, etc., are also capable 
for the excitation of slosh waves in spacecraft fluid systems (Kamotani et al . , 
1981; Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 1993 a,b,c; Hung et al. , 1992 
a,b,c) . 

Among all of the varieties of jitter accelerations listed, accelerations 
induced by slew motion of the spacecraft dominate the forces activated on the 
spacecraft fluid systems. Two coordinate systems (cylindrical and cartesian) 
chosen in this study are (r, $, z) with corresponding velocity components (u r , 
u e> u 2 ) for cylindrical, and (x, y, z) with corresponding velocity components 
(Uj, Uy , u 2 ) for cartesian coordinates. The origin of these two coordinates are 
located at the central bottom of the dewar tank, as shown in Figure 6. The 
spacecraft center of mass, or the geometric center of the spacecraft is located 
at (x c , y c , z c ) — (0, 0, L c ) . The relationships of the coordinate, velocity and 
the force between cartesian and cylindrical coordinates are 


x 


COS0 0 0 


jr 

y 

— 

sin0 0 0 


0 

z 


r-V 

O 

O 


z 


U x 


cos0 -sin0 0 


' u r 

u y 

= 

sin0 cos0 0 


“e 

U z 


0 0 1 
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\Fr) 


cos0 sin© o' 


F x 


= 

-sin0 cos0 0 


F y 

F z\ 


0 0 1 


F ,\ 


(3-29) 


In the derivation of acceleration induced by the slew motion of spacecraft, 
the coordinate system (x", y" , z") is fixed at the spacecraft center of the mass. 
Th© relationships of the coordinate, velocity and acceleration between 
expressions with the origin located at the spacecraft center of the mass (x", y" f 
z") and origin located at the center bottom of the dewar tank (x, y, z) are 


f*xl 
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COS0 0 0 


r 

Ry 

= 

y 

= 

sin0 0 0 


0 

k 


Z ~ L c 


O 

o 
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Z ~ L c 
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' U x 


cos0 -sin0 
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u r 

V y 

= 

u y 

as 

sin0 cos0 
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u e 
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U z 


0 0 
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U z 
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\F X 


COS0 

sin0 

O' 

-1 
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Rz 

= 

F y 

= 

-sin0 

COS0 

0 




slew 

kJ 

l slew 

0 

0 
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F z 

slew 


fk 


cos© sin0 0 


fk 




-sin0 cos0 0 


F y 

(3-33) 

k 

slew 

0 0 1 


F z\ 

slew 


A detailed expression of [R*, Ry, R z ] siew are shown in Equation (3-13) of 
this report. Jitter acceleration is a summation of acceleration induced by slew 
motion and others, such as atmospheric drag on the spacecraft, spacecraft 
attitude motions arising from machinery vibration, thruster firing, crew motion, 
etc. Thus, jitter acceleration can be expressed as 

where f is the jitter frequency (Hz) imposed on the fluid systems of the 
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Ki.rj 

others 
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F z\ 
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others 



cos0 sin0 0 


r f l 
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F r 


= - 

-sin0 cos0 0 


F y 

- 

F 6 

{ 1 •*--isin(27:ft) } (3-34) 


0 0 1 


w 

slew 

F z\ 

2 

others 


spacecraft. 

IV. Non- Inertia Frame Mathematical Formulation of Fundamental Equations 

Dynamical behavior of fluid elements inside the on-orbit spacecraft fluid 
systems are strongly modified by the gravity gradient and gravity jitter 
accelerations. In order to accommodate the impact of gravity gradient 
acceleration, in particular, on the on-orbit fluid motion, one has to consider 
non- inertia frame of the spacecraft bound coordinate rather than adopting inertia 
frame coordinate used in ordinary fluid mechanics formulation. 

Consider a closed circular cylindrical dewar of radius, a, with height, L, 
which is partially filled with cryogenic liquid helium, and rest of the ullage 
is filled with a helium vapor. Angular velocity of rotating cylinder is cj. 
Density and viscosity of liquid helium and helium vapor are p L , /x L , p v , and /Xy, 
respectively (Mason, 1978). Let us use cylindrical coordinates (r, $ t z) , with 
corresponding velocity components (u, v, w) , and corresponding residual gravity 
acceleration, such as gravity gradient components (a gg r , a gg 9 , a gg z ) and gravity 
jitter components (a g j r , a g j a g j z ) . In the derivation of the governing 
equations, accelerations induced by the spinning motion of the spacecraft is 
included in the formulation. The rest of the acceleration such as slew motion, 
atmospheric drag on the spacecraft, spacecraft attitude motions arising from 
machinery vibrations, thruster firing and others, are included in the jitter 
acceleration, shown in Equation (3-34). The governing equations for non- inertia 
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frame of spacecraft bound coordinates can be shown as follows: 
(A) Continuity Equation 


l_d_ 
r dr 


r ad dz 


(4-1) 


(B) Momentum Equations 


/ du du v du v 2 du v do ^ , 

p( lt 17 * 71 e ~~F* , ‘'lz ) af * 2 p “» v * p<a «.'' ,a re.r> 




(4-2) 


/ dv dv v dv uv dv K 1 dp ^ v 

p( lt* u 17*71e*— *“'17 ) "7W -2 p®,u+p - pro. 


♦l*<^-£*£ 3 r 


(4-3) 


„ , dw dw v dw dwx dp , \ ym 

p( 17* u 17*71e*' ,, 17 ) "'17* p<a ^-* a ®- > * p ^‘' <4 ' 4) 


where 


\ 72 - 1 ^ 

r 6r 



+ -i-JL 

r 2 ae 2 


a 2 

dz 2 


(4-5) 


In these formulations, 2cj z v and 2cj z u are the Coriolis acceleration, ra > z 2 is 
the centrifugal acceleration, and ni> z is the angular acceleration induced by the 
spinning motion of the spacecraft. 

In the computation of fluid forces, moment, viscous stress and angular 
momentum acting on the container wall of the spacecraft, one has to consider 
those forces and moment in the inertia frame rather than the non- inertia frame, 
in particular for the case of the spinning motion in the z-axis. To show an 
example, one has to transform those vectors from the non- inertia frame to the 
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inertia frame for the case of spinning motion in the z-axis. 

>'l [coso> z t -sinw z C o] rp ■ 

Fy = sinoo z t coso > z t 0 F y (4-6) 

Fi 0 0 1 . K 

where prime symbol denotes vectors in the inertia frame while those parameters 
without the prime symbol indicate vectors in the non- inertia frame. 

V. Initial and Boundary Conditions of 
Spacecraft Fluid System in Microgravity Environment 
Governing equations of the fluid motion in on-orbit spacecraft fluid 
systems for non- inertia frame of spacecraft bound coordinates have been 
illustrated in Equations (4-1) to (4-6) . These equations shall be combined with 
the characteristics of gravity gradient and gravity jitter accelerations as 
formulated in Equations (3-1) to (3-34). Initial and boundary conditions shall 
be introduced to accommodate solving fluid motion in on-orbit spacecraft fluid 
system for non-inertia frame coordinate (Hung et al. , 1990 a,b,c; 1991 

a,b,c,d,e,f ,g,h, i, j) . 

Let the profile of the interface between gaseous and liquid fluids be given 

by: 

r?(t, r, 0, z) - 0 (5-1) 

The initial condition of the profile of the interface between gaseous and 
liquid fluids at t - t Q is assigned explicitly, and is given by: 

rj(t - t 0 , r, 8, z) - 0 (5-2) 

A set of boundary conditions has to be supplied for solving the equations. 
These initial interface profiles used in this study have been given explicitly 
through the steady state computations made by Hung and Leslie (1988) and Hung et 
al (1989 a,b,c,d) which were checked by the experiments carried out by Leslie 


21 



(1985). These boundary conditions are as follows: 

(1) Along the container wall, the following three boundary conditions 

apply: 

(a) Interface between solid and liquid fluid: No-penetration and no-slip 

conditions assure that both the tangential and the normal components 
of the liquid velocity along the solid walls will vanish. 

(b) Interface between solid and gaseous (vapor) fluid: Similar no- 

penetration and no- slip conditions as that shown for interface between 
solid and liquid fluid will apply. 

-(c) At the location of solid- liquid-gaseous (vapor) three phases 

interface: No-penetration, but not no-slip condition apply. This 

will assure that normal components of liquid and vapor velocities 
along the solid wall vanish, and allow a slipping flow of liquid and 
vapor fluids along the solid wall at three phase interface location. 
The velocity of slipping flow at this location is governed by the 
adhesive forces between fluids (liquid and gaseous) and solid walls. 
Also, at this location of three phase interface, a constant contact 
angle is present in which the behaviors of wet or dry contacts are 
determined by Coulomb interaction between the fluids (liquid and 
vapor) and the surface phenomena (material and roughness) of solid 
walls . 

(2) Along the interface between the liquid and gaseous 

fluids, the following two conditions apply: 

(a) Kinematic surface boundary condition: The liquid (or gaseous) surface 

moves with the liquid (or gas) which implies 
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-Q.O.or 


< 5 - 3 > 

on n ( t=t a , r,B, z) 

(b) Interface stress condition: Across the liquid-vapor interface, the 

stress must be continuous. Based on Landu and Lifshitz (1959), the 
stress across the liquid-vapor interface can be expressed as 

C -(T^) + (5"4) 

where R x and R 2 are the radius of curvatures of two major axes at the point of 
interests on the surface of the liquid-vapor interface. 

The expressions of radius of curvatures R x and R 2 in cylindrical 
coordinates from differential geometry can be shown as 

1 ^ 1 _ EN-2FM+GL _ _ 

X £G-F* (5 ' 5) 

where the relationship of cartesian and cylindrical coordinates for the curved 
surface of liquid-vapor interface is 
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rcos0 

y 

— 

rsin0 

z 


H(t a , r, 0) 


Here, the configuration of the liquid-vapor interface is z - H (t — t a , r, 9 ) ; 



(5-7) 


dxdx dy_ _8y dz dz 

dr ae dr ae dr ae 
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where , 
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Simplifying Equations (5-7) to (5-11), one obtains 


EG-F 2 = r 2 ( 1 + H 2 + -~/f 9 2) = r 2 D 2 


(5-9) 


(5-10) 


(5-11) 


(5-12) 


(5-13) 


(5-14) 
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E = 1 + H 2 , F = H r H 9 , G - I 2 + Hq 2 


IH z ,H 9 ] = [ 


L = 
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{-H a + 


ID 


(r 2 H. 
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-r#ee) 


rD 


and 


, — 

£> = (1 + W r 2 + — flk 2 ) 2 

r 2 

Substituting these relations to Equation (5-4), the radius of the curvature 
on the curved surface of the configuration of liquid-vapor interface in 
cylindrical coordinates can be expressed as follows: 


_1_ 



(r-%) 

r dr D 


— ( — ) ] 
dd id n 


Here, in Equations (5-4) and (5-15) 


(5-15) 


ij 


dut 

dx J 


du J 

dx i 


du ** \ 

SI*"' 


r du k X 


is the viscous stress tensor; /x, the viscous coefficient of the first kind; 
the viscous coefficient of the second kind; P, the pressure; a, the surface 
tension of the liquid- vapor interface; and n j , the unit vector normal to the 
interface; and Sjj, the Kronecker's delta function. Also, subscripts G and L 
denote conditions at gaseous and liquids fluids, respectively, across the liquid- 
vapor interface. 

The fluid stresses across the liquid-vapor interface can be decomposed to 
the components normal (n t , a unit vector) and tangential (t i( a unit vector) to 
the interface. For the component tangential to the interface, one can take a dot 
product of a unit vector tangential to the interface, t if to Equations (5-4) and 


25 



(5-15) , which leads to 


i(T lj t i n J )] L =[(x ij c i n :l )] a ( 5 - 16 ) 

since n t t t - 0. 

For the component normal to the interface, one can also take a dot product 
of a unit vector normal to the interface, n if to Equations (5-4) and (5-15), 
which leads to 


Po-Pl- !<t + 4 <_ ^ )1 (5 * 17) 

For components normal to the interface along the (r, $, z) directions In 
cylindrical coordinates can be obtained by taking dot products of n r , n & , n z 
separately to Equations (5-4) and (5-15), which are expressed as 
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(T ri ^) c - (T XJ nj) L 

< P G ~ P L ) 

n* 

- 

( x 0j n j) g - < T e jOj) L 


n z 
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n z 


( 5 - 18 ) 


where (n r , n dt n 2 ) is the unit vector normal to the interface in cylindrical 
coordinates (r, 9 , z) . 

For a special case of axial symmetry, component normal to the interface, 
shown in Equation (5-17), can be simplified and becomes 

p a - P L - + 


_o d r 

' r dr ' 




(1 + <J> 2 ) 


( 5 - 19 ) 


since d/dd — 0, H r — 3H/3r — dz/dr — H$ — 0 and D — (1 + ^ 2 ) 1/2 for the case 
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of axial symmetry. 

VI . Characteristics of Slosh Wave Induced Fluctuations 
in Fluid Moment and Angular Momentum 

Slosh wave induced fluctuations in the fluid system of the rotating dewar 
introduce time -dependent disturbances in moment and angular momentum of 
spacecraft fluid system. In this study, there are induced angular velocities 
along the yawing, pitching and rolling axes due to the fluid motion inside 
rotating container. These angular velocities in yawing, pitching and rolling 
axes, caused by the fluid flows in a partially liquid- filled container, readjust 
the angular velocity in rolling axis. 

In order to accommodate the spacecraft dynamics of yawing, pitching and 
rolling, cylindrical coordinates (shown in Figure 2) of rotating container is 
transformed into cartesian coordinates based on (x, y, z) — (rcos0, rsinifl, z) 
with corresponding velocity components (V x , V y , V z ) - (ucosfl - vsinfl, usinfl + 
vcosfl, w) . If spacecraft is rotated with respect to mass center at (r c , 0 C , z c ) 
in cylindrical coordinates, location of mass center in cartesian coordinates 
becomes (x c , y c , z c ) * (r c cos0 CJ r c sin0 c , z c ) . Induced angular velocities (w x , a> y , 
w z ) in cartesian coordinates becomes 


G) * 

X 


K — K K 

xx *^xy 



G> 

y 

= 

K K — ir 

1K yx ^yy , yz 


v r 



-K zx K zy K zz 


[vj 


where 


K^K^K^O 


27 



xy 


y* 


(z-z c ) 


[ (y-y c ) 2 + (z-z c ) 2 ] _1 


[ (x-x c ) 2 + (z-z c ) 2 ] _1 


K x 

K ' 


= (y-y e ) 


[ (y-y c ) 2+ (z-z c ) 2 ] _1 
[ (x-x c ) 2 +(y-y c ) 2 ] _1 


yz 


zy. 


= (x-xJ 


[ ( x-x c ) 2 + (z-z c ) 2 ] _1 
[ (x-x c ) 2 + (y-y c ) 2 ] ' 1 


As the velocity components are given by 
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the relationship between the components of induced angular velocity and flow 
velocity in cylindrical coordinates can be expressed in the following 
formulation: 
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The moment of spacecraft can be computed from the time rate of change of 
the angular momentum, i.e., 
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(6-5) 
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where - (w x , <j y , u> 2 ) denotes pitching , yawing and rolling angular velocities 
of spacecraft in inertia frame. 

VII. Mathematical Formulation of Fluid Stresses 
and Moment Fluctuations Due to Slosh Waves 
For the purpose of considering large amplitude slosh wave activated fluid 
stresses exerted on the solid walls of the dewar, the fluid stresses are 
decomposed into the tangential and normal components to the walls which can be 
expresses as follows: 


tt ( du_ du fl \ ~ ^ 


(7-1) 


where II t denotes the tangential component of fluid stresses; 1^, the normal 
component of fluid stresses; P, the thermodynamic pressure; u^, fluid velocity 
in a direction; t a , unit vector tangential to the wall; n^ f unit vector normal 
to the wall; the molecular viscosity coefficient of fluid; and 5^, the 
Kronecker's delta function. Subscripts a and 0 imply the directions of flow 
fields . 

Figures 6(A) and 6(B) show the geometry of GP-B dewar propellant tank in 
both r-z and r-0 planes, respectively. In order to make the computation of fluid 
stresses match the geometry of the dewar tank, mathematical formulations have 
been divided into three sections: (A) Top wall (dome) section, (B) Bottom wall 


dXa dx t 


* A A 


(7-2) 
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(dome) section, (C) Cylindrical section (including probe section of the inner 
wall of dewar) and (D) Baffle plates section. There are several plates of baffle 
inserted in the dewar [see Figure 7(A) and 7(B) for baffle board illustration]. 
Baffles with the shape of hollow circular plate with an inner radius R lf and an 
outer radius R 2 are installed along the probe column of the dewar located at z 
- L if where i - 1, 2, n plates and the thickness of each plate is d. 

(A) Top Wall (Dome) Section: 



(n c ) Top + "§r) cos ^ 

(7-3) 


(7-4) 



(7-5) 

(B) Bottom Wall (Dome) 

Section: 



<n c > a. h-4?w°s2* 


dz dr / 


TT J 0 

Bottom 
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l(T^ + a?) cos ^lTle + ai) sin<,> J 


(7-6) 
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Bottom Wall ^ ^) sin2 * < 7 - s > 

(C) Cylindrical Section (including probe section of the inner wall of dewar): 


(n c )£fi* dricaJ =^||+|^) ( 7 - 9 ) 

(n c )^ indricaJ =n(^|H + |x| (7-io) 
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(7-11) 


<n„> 
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(D) Baffle Plates Section: 


(II ) r ' z = U (— + 4^) 

c Baffle ** dz dl 


(7-12) 
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(7-14) 


where <£ is the azimuth angle of the dome; and (II t ) r z and (II t ) r 9 denote tangential 
stresses in r-z and r-$ planes of the dewar, respectively. Velocity components 
in cylindrical coordinates of (r, 9, z) are shown as (u, v, w) . 

The stress distribution shown in Eqs . (7-3) to (7-14) can be integrated 
with respect to area and obtain the tangential and normal forces on top wall, 
bottom wall, and cylindrical sections of the dewar. 

(A) Top Wall (Dome) Sections: 
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(B) Bottom Wall (Dome) Section: 
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(C) Cylindrical Section (including probe section of the inner wall of dewar) 
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(D) Baffle Plate Section: 
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Baffle 


where and <f> denote the radius and the azimuth angle of the dome section, 
respectively; and R cl and R c2 are the radii of the inner and outer walls of the 
cylindrical section of the dewar (see Figs. 6 and 7). 

In order to accommodate the spacecraft dynamics of pitching, yawing and 
rolling, cylindrical coordinates (shown in Figs. 6 and 7) of the rotating 
container is transformed into cartesian coordinates based on (x, y, z) - (rcos# , 
rsintf, z) with corresponding velocity components (v x , v y , v z ) - (ucosfl - vsin0, 
usin0 + vcos0, w) . For the case of the AXAF-S spacecraft, the axes of slew 
motion are always fixed at the mass center of the spacecraft which is located at 
(x c , y c , z c ) - (0, 0, L c ) where L c is the height of the axis of slew motion (see 
Figure 4). To fulfill this goal, stress distributions, shown in Eqs . (7-15) to 
(7-18), have to be recalculated in the (x, y, z) directions (see Figs. 6 and 7). 
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'( [ //“ (IIc) SS*rcyJiud.rSin0 *//< n J InnerCyllnder < ~ O S 0] i? cl d8dz} 
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^OB^3in^dlpde}^Jf (n t ) I I ^ mrC yllJ > dmr^cJ dedz * [[(^c) r oHmrCynndmT^c2^ dz ] (7“2l) 

The moment of stress force acting on the dewar wall of the container can 
be computed from the cross product of moment arm, [which is the perpendicular 
distance from the location of proof mass of the spacecraft to the total forces] , 
and the total forces. Components of the moment can be computed from the 
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following formulations: 
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+ [f / < n J top waii + <n fl ) SoCCOffl „ aii ]i?dCos 3 <J)-cos0-sined(t>d9 
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where denotes the height of the dome (see Fig. 2). 

After the integration of Eqs . (7-19) thru (7-24), components of the moment 
can be rewritten as the following expression: 
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where L*, ly , and denote the components of the moment arm along x, y, and z 
axes, respectively. 

By using the relation of F a L* - 0, or F x L* + F y Ly + F z (Lj - L c ) - 0, 
moment arms of the moment of fluid stress moment induced by the slosh wave 
excitation can be computed from the following relations: 

L x ] \FJi z -F& 

L y = 2 h i F z M x~ F x M z (7-26) 

L z -L c \ F*+F*+F* [F^-FyM x \ 

VIII. Methods of Numerical Simulation 
Detailed descriptions of the computational algorithm applicable to 
cryogenic fluid management under microgravity are also given in our earlier 
studies (Hung et al . , 1990 a,b,c). In this report, a full-scale AXAF-S 

spacecraft propellant dewar tanks with a radius of 68 cm and a height of 145 cm 
will be used in the numerical simulation. The propellant tank is 70% filled with 
cryogenic liquid helium and the rest of the ullage is filled with helium vapor. 
The temperature of cryogenic helium is 1.3 K. In this study the following data 
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were used: liquid helium density — 0.145 g/cm 3 , helium vapor density - 0.00147 
g/cm 3 , fluid pressure - 1.66 x 10 3 dyne/cm 2 , surface tension coefficient at the 
interface between liquid helium and helium vapor — 0.346 dyne/cm, liquid helium 
viscosity coefficient - 1.12 x 10‘ 4 cm 2 /s ; and contact angle - 5°. The initial 
profiles of the liquid-vapor interface for the rotating dewar are determined from 
computations based on algorithms developed for the steady state formulation of 
microgravity fluid management (Hung et al . , 1990 a,b,c). 

A staggered grid for the velocity components is used in this computer 
program. The method was developed by Harlow and Welch (1965) for their MAC 
(marker-and-cell) method of studying fluid flows along a free surface. The 
finite difference method employed in this numerical study was the "Hybrid Scheme" 
developed by Spalding (1972) . The formulation for this method is valid for any 
arbitrary interface location between the grid points and is not limited to middle 
point interfaces (Patankar and Spalding, 1972). An algorithm for a semi- implicit 
method (Patankar, 1980) was used as the procedure for modeling the flow field. 
The time step is determined automatically based on the size of the grid points 
and the velocity of flow fields. A detailed description of the computational 
algorithm applicable to microgravity fluid management is illustrated in our 
earlier studies (Hung et al . , 1990 a,b,c). Figures 8(A) and 8(B) show the 
similar distribution of grid points for the dewar tank for the AXAF-S spacecraft 
in the radial -axial plane and radial -circumferential plane, respectively, in 
cylindrical coordinates. 

IX. Spacecraft Sloshing Dynamics Associated With Spinning and/or Slew Motions 

By using the mathematical formulations illustrated in Sections III and IV, 
subject to the initial and boundary conditions in Section V, spacecraft sloshing 
dynamics associated with slew motions depending upon the specific scientific 
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missions assigned to the spacecraft are investigated. Two examples are given to 
illustrate sloshing dynamics driven by each of the gravity gradient and jitter 
accelerations associated with slew motion for the AXAF-S spacecraft. 

(IX-A) Spacecraft Sloshing Dynamics Driven by Gravity Gradient Acceleration 
Associated With Slew Motion 

As we indicated earlier, AXAF-S spacecraft is capable to observe point and 
extended sources of active galactic nuclei, clusters of galaxies, supernova 
remnants, x-ray binaries, etc., through spacecraft slew motion of pointing 
control. Assume that slew motion is along the y"-axis (see Figure 1), gravity 
gradient acceleration associated with slew motion can be computed from the non- 
inertia frame expressions of Equations (3-1) and (3-26). It is assumed that the 
slew motion operates in 90° in 10 minutes. Non- inertia frame fundamental 
equations, shown in Equations (4-1) to (4-5) shall be adopted with making u z and 
o) z equal to zero because there is no spinning motion in the z-axis. Initial and 
boundary conditions, shown in Equations (5-1) to (5-19), shall also be considered 
in our computation. 

In this example, spacecraft sloshing dynamics driven by gravity gradient 
acceleration associated with slew motion in the y"-axis, shown in Figure 1, have 
been investigated. As the orbital period of AXAF-S spacecraft is 97.6 min and 
period of slew motion in the y ,f -axis is 600 s, the component of gravity gradient 
acceleration along the (x, y, z) directions acted on the fluid mass located at 
(r, 6, z) - (12 cm, n/2 , 3 cm) is shown in Figure 9. This figure shows that the 
magnitude of gravity gradient acceleration is on the order of 1CT 7 ^ for AXAF-S 
dewar on its operation orbit. The distance from the spacecraft mass center to 
the bottom of the dewar, L c , shown in Figure 1, is 257.8 cm. The grid point 
generation of the dewar tank, shown in Figure 8, is adopted in this example of 
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the AXAF - S dewar . 


Gravity gradient acceleration indicates that acceleration acted on any 
fluid mass inside the container increases two units of acceleration per unit of 
distance measured from the spacecraft mass center to the location of the fluid 
mass parallel along the radial axis from the spacecraft mass center to the center 
of the Earth (parallel to unit vector r c shown in Figure 1) while the 
acceleration acted on the fluid mass decreases one unit of acceleration per unit 
of the shortest distance measured from the location of the fluid mass to the 
radial axis along the vector from the mass center of the spacecraft to the center 
of the Earth (Forward, 1982) . As the magnitude and the direction of gravity 
gradient acceleration acted on each fluid mass is strongly dependent upon how far 
the location of the fluid mass deviates from the spacecraft mass center measured 
along the axis parallel to the vector r c which varies with respect to time, it 
can be anticipated that gravity gradient acceleration acted on the fluid mass is 
different for fluid mass at different locations in the container. 

The equilibrium shape of the liquid-vapor interface for a dewar with 70% 
liquid- filled level under a residual gravity environment below 10‘ 7 g^ as that 
shown in Figure 9, is a sphere. Figure 10(A) shows the initial shape of the 
interface in the r-z plane at $ - 0° and 180°; Figure 10(B) shows the initial 
profile of the liquid-vapor interface in the r-z plane at 6 - 90° and 270°; 
Figure 10(C) shows the initial profile of the liquid-vapor interface in the r-0 
plane at height z - 95.9 cm; and Figure 10(D) shows the initial profile of three- 
dimensional liquid-vapor interface . 

Figure 11 shows the time sequence evolution of the three-dimensional 
dynamical behavior of the interface oscillations driven by gravity gradient 
acceleration associated with slew motion. For the convenience of comparison. 
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figures of liquid-vapor interface profiles with the same values of the time 
sequences chosen throughout this paper are at time t - 334, 392, 431, 456, 487, 
524, 554, 588, 600, 695, 784 and 800 seconds. It clearly shows that there are 
a series of asymmetric oscillations excited along the surface of sloshing 
dynamics governed liquid-vapor interface driven by asymmetric gravity gradient 
acceleration associated with slew motion. 

The evolution of the sloshing dynamics governed interface oscillations at 
various cross-sections driven by gravity gradient acceleration associated with 
slew motion are also examined. Figures 12 and 13 show the time sequence of the 
sloshing dynamics governed liquid-vapor interface profiles, driven by the same 
gravity gradient acceleration associated with slew motion as that shown in Figure 
11, in the vertical r-z plane at 9 - 0° and 180°, and at 6 - 90° and 270°, 
respectively. It indicates that the spherical shape bubble (helium vapor) 
configurations change from axial symmetric to asymmetric profiles in both planes 
during the course of slew motion. 

Figure 9 shows that gravity gradient acceleration associated with slew 
motion is pointing toward the southwestern direction in the x-z plane (or r-z 
plane at 6 - 0° and 180°) , and also is pointing toward the southern direction in 
the y-z plane (or r-z plane at 9 - 90° and 270°). These indicate that liquid is 
pushed toward the southwestern direction (the bubble is pushed toward the 
northeastern direction) in r-z plane at 8 - 0° and 180° while the liquid is 
pushed toward the southern direction (the bubble is pushed in the northern 
direction) in r-z plane at 9 - 90° and 270°. Figures 12 and 13 exactly 
demonstrate these results. 

Figure 14 shows the time evolution of the sloshing dynamics governed 
liquid-vapor interface oscillations driven by the same gravity gradient 
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acceleration associated with slew motion, as that shown in Figure 11, in the 
horizontal (r-0) plane for height at z — 95.9 cm. It shows that a circular- 
shaped bubble configuration of the liquid-vapor interface changes from axial 
symmetry [see Figure 10(C)] to asymmetry (see Figure 14) profiles. Based on 
Figure 9, it shows that gravity gradient acceleration associated with slew motion 
is pointing toward the western direction in the x-y plane (or r-0 plane). This 
implies that the liquid is pushed toward the western direction while bubble is 
pushed toward the eastern direction in the r-0 plane. In the meanwhile, the 
bubble is also pushed in a upward direction along the spinning axis of the dewar 
as that shown in Figures 12 and 13. This scenario exactly reflects the evolution 
of the bubble shown in Figure 14 that the bubble is pushed toward the eastern 
direction and also gradually becomes smaller due to the upward movement of the 
bubble . 

Comparison of Figures 9,11,12,13 and 14 illustrate some peculiar behavior 
of cryogenic helium II fluids with temperature below the A -point (2.17 K) in 
which liquid helium demonstrates a number of remarkable properties of 
superfluidity such as extremely low viscous and surface tension coefficients 
reacted to the disturbances driven by gravity gradient acceleration associated 
with slew motion. These results can be concluded as follows: (a) Both negative 
direction of gravity gradient acceleration in x and z components imply that the 
liquid is pushed toward negative x and z directions while bubble is pushed toward 
positive x and z directions. (b) Comparison of the magnitudes of gravity 
gradient acceleration in the x, y and z components, shown in Figure 9, indicated 
that g 2 > g x > g y . This magnitude of driving force reflects the results of the 
bubble mass center fluctuations, shown in Figures 11, 12 and 13. (c) 
Configuration of bubble fluctuations at time t - 334 s, shown in Figures 11 and 
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12, indicate a near perfect shape of a sphere, as that shown in Figure 10, even 
though an extensive gravity gradient acceleration has continuously applied to the 
fluid system starting at t - 0 , shown in Figure 9. (d) Uneven and imbalanced 
flow velocities toward the southwestern direction create similar uneven and 
imbalanced pressure distribution reacted to the bubble toward the northeastern 
direction. Because of extremely low surface tension coefficients of helium II 
between the liquid-vapor interface, a deformed irregular concave and convex- 
shaped oscillating bubble is created. 

Figure 15 shows time fluctuations of the locations of bubble mass centers 
of the fluids inside the dewar container due to sloshing dynamics driven by 
gravity gradient acceleration associated with slew motion. The values of bubble 
mass center fluctuations are (Ax c , Ay c , Az c ) - (11.8, 1.25, 18.4) cm. It shows 
Az c > Ax c > Ay c for bubble mass center fluctuations driven by gravity gradient 
acceleration associated with slew motion. 

The trend of the bubble center fluctuations, Az c > Ax c > Ay c , exactly 
reflect the values of major driving forces of gravity gradient acceleration 
associated with slew motion in which it shows g 2 > g x > g y as that shown in 
Figure 9 . 

(IX-B) Spacecraft Sloshing Dynamics Driven by Jitter Acceleration Associated 
With Slew Motion 

By using the mathematical formulations illustrated in Sections III to IV, 
one can numerically simulate spacecraft sloshing dynamics associated with 
spinning and/or slew motions depending upon the specific scientific missions 
assigned to the spacecraft. As we indicated earlier, AXAF-S spacecraft is 
capable to observe point and extended source of active galactic nuclei, clusters 
of galaxies, supernova remnants, x-ray binaries, etc., through spacecraft slew 
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motion of pointing control. 

If slew motion operates with a range of 90° in 10 minutes (- 600 s) , the 
component of jitter acceleration, based on Equations (3-1) to (3-34), along the 
(x, y, z) directions acted on the fluid mass located at (r, 9, z) - (12 cm, n/2, 

3 cm) is shown in Figure 16. This figure shows that the magnitude of gravity 

gradient acceleration is on the order of 10 5 go- 

The equilibrium shape of the liquid-vapor interface for a dewar with 70% 
liquid- filled level under a residual gravity environment below 10' 6 go is a 
sphere. Figure 10(A) shows the initial shape of the interface in the r-z plane 
at 9 - 0° and 180°; Figure 10(B) shows the initial profile of the liquid-vapor 
interface in the r-z plane at 9 - 90° and 270°; Figure 10(C) shows the initial 
profile of the liquid-vapor interface in the r -9 plane at height z - 95.9 cm; and 
Figure 10(D) shows the initial profile of three-dimensional liquid-vapor 
interface . 

Figure 17 shows the time sequence evolution of the three-dimensional 
behavior of the interface oscillations driven by jitter acceleration associated 
with slew motion. It is shown in this figure that a time sequence evolution of 
liquid-vapor interface profiles at time t - 60.4, 258, 323, 354, 379, 403, 430, 
528, 572, 628, 663, and 800 seconds are illustrated. It clearly shows that there 
are a series of asymmetric oscillations excited along the surface of sloshing 
dynamics governed liquid-vapor interface driven by asymmetric jitter acceleration 
associated with slew motion. 

Figure 18 shows the time sequence of the sloshing dynamics governed liquid- 
vapor interface profiles, driven by the same jitter acceleration associated with 
slew motion as that shown in Figure 17, in the vertical r-z plane at 9 - 0° and 
180°. It clearly indicates that the spherical shape bubble (helium vapor) 
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configurations change from axial symmetric to asymmetric profiles during the 
course of slew motion. 

The evolution of the sloshing dynamics governed interface oscillations at 
various cross-sections driven by gravity jitter acceleration associated with slew 
motion are also examined. Figures 18 and 19 show the time sequence of the 
sloshing dynamics governed liquid-vapor interface profiles, driven by the same 
gravity jitter acceleration associated with slew motion as that shown in Figure 
17, in the vertical r-z plane at 6 - 0° and 180°, and at $ - 90° and 270°, 
respectively. It indicates that the spherical shape bubble (helium vapor) 
configurations change from axial symmetric to asymmetric profiles in both planes 
during the course of slew motion. 

Figure 16 shows that gravity jitter acceleration associated with slew 
motion is pointing toward the southeastern direction at the very beginning and 
then is pointing toward the southwestern direction at the end in the x-z plane 
(or r-z plane at 6 - 0° and 180°) , and also is pointing toward the southern 
direction in the y-z plane (or r-z plane at ^ - 90° and 270°) . These indicate 
that liquid is pushed toward the southeastern and then toward the southwestern 
directions (the bubble is pushed toward the northwestern and then toward the 
northeastern directions) in r-z plane at 9 - 0° and 180° while the liquid is 
pushed toward the southern direction (the bubble is pushed toward the northern 
direction) in r-z plane at J - 90° and 270°. Figures 18 and 19 exactly 
demonstrate these results. 

Figure 20 shows the time evolution of the sloshing dynamics governed 
liquid-vapor interface oscillations driven by the same gravity jitter 
acceleration associated with slew motion, as that shown in Figure 17, in the 
horizontal (r -0) plane for height at z — 95.9 cm. It shows that a circular- 
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shaped bubble configuration of the liquid-vapor interface changes from axial 
symmetry [see Figure 10(C)] to asymmetry (see Figure 20) profiles. Based on 
Figure 16, it shows that gravity jitter acceleration associated with slew motion 
is pointing toward the eastern direction at the beginning and then pointing 
toward the western direction at the end in the x-y plane (or r-0 plane). This 
implies that the liquid is pushed toward the eastern direction and then toward 
the western direction while bubble is pushed toward the western direction and 
then toward the eastern direction in the r-0 plane. In the meanwhile, the bubble 
is also pushed toward the upward direction along the spinning axis of the dewar 
as that shown in Figures 18 and 19. This scenario exactly reflects the evolution 
of the bubble shown in Figure 20. 

Figure 21 shows time evolution of the locations of bubble mass centers of 
fluids inside the dewar container due to sloshing dynamics driven by jitter 
acceleration associated with slew motion. The values of bubble mass center 
fluctuations are (Ax c , Ay c , Az c ) - (28.9, 0.44, 30.2) cm. It shows Az c > Ax c > 
Ay c for bubble mass center fluctuations driven by jitter acceleration associated 
with slew motion. Behavior of bubble mass center fluctuations, shown in Figure 
21 are the exact reflection of the behavior of jitter accelerations, shown in 
Figure 16. 

Comparison between Figures 16 to 21 illustrate some peculiar behavior of 
cryogenic helium fluids with temperature below A-point (2.17 K) in which helium 
demonstrates a number of remarkable properties of superfluidity such as extremely 
low viscous and surface tension coefficients reacted to the disturbances driven 
by jitter acceleration associated with slew motion. It can be concluded as 
follows: (a) Jitter acceleration associated with slew motion started at time t 
— 0. However, bubble mass center fluctuations did not start to react to the 
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driven force until t - 80 s for x c and t - 160 s for z c . In other words, bubble 
configuration was in perfect spherical-shaped at 80 s after jitter acceleration 
associated with slew motion was applied. (b) x-component jitter acceleration 
associated with slew motion was applied to the fluid element with positive value 
from t - 0 to 300 s; negative value from t - 300 to 600 s; and zero value after 
t - 600 s. However, westward movement of bubble continued from t - 80 to 450 s, 
and then switched to eastward movement to t - 800 s which was 200 s after 
acceleration vanished. Obviously, there is a phase-shift between action of force 
and reaction of motion, (c) The z-component jitter acceleration associated with 
slew motion was applied to the fluid element with negative value from t — 80 to 
500 s and the zero value after t - 500 s. However, the northward movement of the 
bubble started at t - 160 s; bounced back from the wall at t - 530 s; and then 
the southward movement continued to t - 800 s which was 300 s after acceleration 
vanished. Obviously, the motion continued for a long period of time even after 
the applied force vanished due to extremely low viscosity of helium II fluids, 
(d) An intensive oscillation of bubble with a deformation of irregular concave 
and convex- shaped continued and sustained for several hundred seconds after the 
applied force vanished due to extremely low surface tension coefficient for 
helium II fluids. 

X. Sloshing Dynamics Induced Angular Momentum and Moment Fluctuations 
Driven by Gravity Gradient and Jitter Accelerations 
Associated With Slew Motion 

(X-A) Angular Momentum and Moment Fluctuations Driven by Gravity Gradient 

Acceleration Associated With Slew Motion 
Figure 22 shows the computed time variation of the fluctuations of angular 
momentum driven by gravity gradient acceleration associated with slew motion. 
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This figure shows the following results: (a) The values of angular momentum 
fluctuations are (AH X , AHy, AH Z ) - (0.032, 2.602, 0.015) 10 6 gcm 2 /s, it clearly 
indicates AHy > AH X > AH 2 . The maximum absolute values of angular moment are Max 
(|Hj, |Hy|, |Hj) - (0.032, 2.602, 0.015) 10 6 gcn 2 /s. It also indicates | Hy | 
> |Hj > | H z | . (b) The initial values of Hj, Hy and H x start from zero value in 
non- inertia frame. (c) Variations of Hy is very much following the trend of 
angular displacement in the y-axis due to the slew motion, shown in Figure 9, 
while the variations of H* and H 2 are much smaller than Hy because there is no 
slew motion applied on the x- and z-axes. 

Figures 23(A), 23(B) and 23(C) show variations of fluid moments due to 
sloshing dynamics driven by gravity gradient acceleration associated with slew 
motion along the x, y and z axes, respectively. The values of fluid moment 
fluctuations are (AM^, AMy, AM Z ) - (613.1, 8235.4, 129.0) dyne cm. The maximum 
absolute values of fluid moment are Max (Im*!, | My | , |M Z |) - (412.4, 6101.1, 
69.2) dyne* cm. It shows AMy > AM X > AM Z and l«yl > |Mxl > |M 2 | . Characteristics 
of the fluctuations of fluid moments due to sloshing dynamics driven by gravity 
gradient acceleration associated with slew motion in the y-axis draw three point 
conclusions similar to that drawn for the fluctuations of angular momentum of 
fluid system except that the trend of fluid moment fluctuations follow that of 
angular velocity of the spacecraft slew motion in the y-axis, shown in Figure 4, 
rather than that of angular displacement for angular momentum because fluid 
moment is time derivative of angular momentum while the angular velocity is the 
time derivative of angular displacement. 

(X-B) Angular Momentum and Moment Fluctuations Driven by Gravity Jitter 
Acceleration Associated With Slew Motion 

Figures 24(A, 24(B) and 24(C) show the computed time variation of the 
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fluctuations of angular momentum driven by gravity jitter acceleration associated 
with slew motion along the x, y and z axes, respectively. This figure shows the 
following results: (a) The values of angular momentum fluctuations are (AH,, 
AHy, AH,) - (0.4, 17.55, 0.161) 10 6 gxm 2 /s, it clearly indicates AHy > AH, > AH Z . 
The maximum absolute values of angular moment are Max ( I H* I > I Hy I • I H z I ) " 
(0.275, 8.92, 0.091) 10 6 g cm 2 /s. It also indicates |Hy| > |H*| > |hJ . (b) The 
initial values of H, , Hy and H x start from zero value in non- inertia frame. (c) 
Variations of Hy dominates the whole spectrum of angular momentum fluctuations 
because angular displacement in the y-axis is the only slew motion in the entire 
process, while the variations of H, and H z are much smaller than Hy because there 
is no slew motion applied on the x- and z-axes. 

Figures 25(A), 25(B) and 25(C) show variations of fluid moments due to 
sloshing dynamics driven by gravity jitter acceleration associated with slew 
motion along the x, y and z axes, respectively. The values of fluid moment 
fluctuations are (AM,, AMy, AM Z ) - (13.38, 200.79, 4.32) 10 3 dynecm. The maximum 
absolute values of fluid moment are Max (|M,|, | My | , | M, | ) - (6.77, 100.75, 2.33) 
10 3 dyne cm. It shows AMy > AM, > AM Z and | My | > |M,| > | M z | . Characteristics 
of the fluctuations of fluid moments due to sloshing dynamics driven by gravity 
jitter acceleration associated with slew motion in the y-axis draw three point 
conclusions similar to that drawn for the fluctuations of angular momentum of 
fluid system. 

XI. Characteristics of Sloshing Dynamics Induced Fluid Stresses 
and Moment Fluctuations Exerted on the Dewar Container 

In this study, characteristics of the asymmetric fluctuations of liquid- 
vapor interface due to sloshing dynamics driven by each of the gravity gradient 
and jitter accelerations associated with slew motion were discussed in Section 
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IX* In Section VIII, the mathematical formulation of the fluctuations of fluid 
stresses and the moment of fluid stresses exerted on the dewar container have 
been illustrated. With references to Sections VIII and IX, one can calculate the 
fluid stresses and the associated moments exerted on the dewar container of 
spacecraft dewar tank. 

(XI-A) Fluid Stress Forces and Moment Fluctuations Exerted on Dewar Container 
Driven by Gravity Gradient Acceleration Associated With Slew Motion 

Figure 26(A), 26(B) and 26(C) show the computed time variation of the 
fluctuations of fluid stress forces exerted on the dewar container driven by 
gravity gradient acceleration associated with slew motion in the y-axis with 
components along the x- , y- , and z-axes, respectively. These figures show the 
following results: (a) The values of fluid stress force fluctuations are (AF X , 
AF y , AF Z ) - (59.2, 4.11, 61.02) dynes, it clearly indicates AF Z > AF X > AF y . The 
maximum absolute values of fluid stress force are Max (|F X |, |F y |, |F 2 [) - 
(47.42, 4.11, 78.68) dynes. It also indicates |F 2 | > |F X | > |F y |. (b) The 
initial values of F x and F y start from zero value while that of the F z starts 
from non-zero value. This is due to the fact that the major driving force of 
gravity gradient acceleration associated with slew motion at the beginning of 
slew motion is equal to zero along the x- and y-axes, and is non-zero along the 
z-axis. (c) The characteristics of fluid stress forces and their fluctuations 
are more likely following the trend of the major driving force of gravity 
gradient acceleration associated with slew motion in the y-axis as that shown in 
Figure 9. Comparison between Figures 26 and 9 reflect that the fluid system 
serves as a damping modulator for the forces acting on the fluid flows and then 
responds back to the spacecraft. 

Figures 27(A), 27(B) and 27(C) show the variations of fluid stress moments 
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exerted on the dewar container due to sloshing dynamics driven by gravity 
gradient acceleration associated with slew motion with components along the x, 
y and z axes, respectively. These figures show the following results: (a) The 
values of fluid stress moment fluctuations are (AM*, AMy, AM*) - (770.6, 10784.5, 
0.0004) dyne cm. The maximum absolute values of fluid stress moment are Max 
(|M*|, | My | , |M Z |) - (770.6, 8595.9, 0.0004) dyne*cm. (b) As slew motion is in 
the y-axis, this makes the magnitudes of both AMy and Kl to be the maximum, 
(c) As slew motion is in the y-axis, this makes the value of moment arm along the 
y-axis to be zero; and also as the y-component of the major driving force of 
gravity gradient acceleration associated with slew motion in the y-axis is equal 
to zero, this makes the magnitude of |F y | to be near zero. Since both | F y | — 
Ly - 0, it leads M z - L*F y - LyF x - 0. 

Figure 28(A), 28(B) and 28(C) show time fluctuations of moment arms of 
fluid stress moments exerted on the dewar container due to sloshing dynamics 
driven by gravity gradient acceleration associated with slew motion in the y- 
axis for components along the x- , y- , and z-axes, respectively. These figures 
show the following results: (a) The values of moment arm fluctuations are (AL^, 
ALy, AL*) - (119.2, 11.7, 156.1) cm. The maximum absolute values of moment arm 
are Max ( | L* | , | Ly | , | | > — (91 . 3 , 11.7, 257 .8) cm. It shows AL* > AL* > ALy and 
|l*| > |Lj > IM- (b) As spacecraft slew motion is along the y-axis, this 
makes the values of Ly to be the minimum, and My to be the maximum. (c) As L z 
is the function of F^My, the values of My is the maximum and F x is near to the 
maximum. This makes the values of L z to be the maximum. (d) The trend of the 
fluctuations of momentum arm is very much like that of the fluid stress forces 
acting on the dewar container of the spacecraft. Again, the trend of the 
fluctuations of moment arm is very much like that of the major driving forces of 
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gravity gradient acceleration associated with slew motion in the y-axis (see 
Figures 28 and 9). 

(XI-B) Fluid Stress Forces and Moment Fluctuations Exerted on Dewar Container 
Driven by Gravity Jitter Acceleration Associated With Slew Motion 

Figure 29(A), 29(B) and 29(C) show the computed time variation of the 
fluctuations of fluid stress forces exerted on the dewar container driven by 
gravity jitter acceleration associated with slew motion in the y-axis with 
components along the x- , y- , and z-axes, respectively. These figures show the 
following results: (a) The values of fluid stress force fluctuations are (AF X , 
AF y , AF Z ) - (2.27, 0.04, 1.37) 10 3 dynes, it clearly indicates AF X > AF 2 > AF y . 
The maximum absolute values of fluid stress force are Max (|F X |, |F y |, |F Z |) - 
(1.15, 0.02, 1.20) 10 3 dynes. It also indicates |F 2 | > |F X | > |F y |. (b) The 
initial values of F x , F y and F z all start from zero value because there was no 
slew motion applying to the system at time t - 0. (c) Magnitudes and 
fluctuations of F y is much smaller than F x and F z because the major driving 
moment of slew motion is along the y-axis. (d) The characteristics of fluid 
stress forces and their fluctuations are more likely following the trend of the 
major driving force of gravity jitter acceleration associated with slew motion 
in the y-axis as that shown in Figure 16. Comparison between Figures 29 and 16 
reflect that the fluid system serves as a damping modulator for the forces acting 
on the fluid flows and then responds back to the spacecraft. The small 
fluctuations shown in the forces response back to the dewar container are due to 
the disturbances caused by the sloshing dynamics of the fluid system. 

Figures 30(A), 30(B) and 30(C) show the variations of fluid stress moments 
exerted on the dewar container due to sloshing dynamics driven by gravity jitter 
acceleration associated with slew motion with components along the x, y and z 
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axes, respectively. These figures show the following results: (a) The values 

of fluid stress moment fluctuations are (AM*, AMy, AM X ) - (7.53, 422.79, 0.0001) 
10 3 dynecm. The maximum absolute values of fluid stress moment are Max 
| My | , |M X | ) - (4.31, 214.53, 0.0006) 10 3 dyne cm. (b) As slew motion is in the 
y-axis, this makes the magnitudes of both AMy and to be the maximum. (c) 

As slew motion is in the y-axis, this makes the value of moment arm along the y- 
axis to be zero; and also as the y-component of the major driving force of 
gravity jitter acceleration associated with slew motion in the y-axis is equal 
to zero, this makes the magnitude of | F y | to be near zero. Since both |F y | ~ 

Ly - 0, it leads M* - LyF z - L^Fy - 0 and M z - L*Fy - LyF x - 0. 

Figure 31(A), 31(B) and 31(C) show time fluctuations of moment arms of 
fluid stress moments exerted on the dewar container due to sloshing dynamics 
driven by gravity jitter acceleration associated with slew motion in the y-axis 
for components along the x- , y- , and z-axes, respectively. These figures show 
the following results: (a) The values of moment arm fluctuations are (AL*, ALy, 

AL*) - (190.4, 56.1, 186.6) cm. The maximum absolute values of moment arm are 
Max (|l*i. mi. mi) - (95.4, 49.8, 186.6) cm. It shows AL* > AL* > ALy and |Lj| 
> mi > mi- (b) As spacecraft slew motion is along the y-axis, this makes the 
values of Ly to be the minimum, and My to be the maximum. (c) As 1^ is the 

function of F*My (FyM* - 0 because F y - 0) , the values of My is the maximum and 

F x is near to the maximum. This makes the values of L* to be the maximum. (d) 
The trend of the fluctuations of momentum arm is very much like that of the fluid 
stress forces acting on the dewar container of the spacecraft. Again, the trend 
of the fluctuations of moment arm is very much like that of the major driving 
forces of gravity jitter acceleration associated with slew motion in the y-axis 
(see Figures 31 and 16). 
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XII. Discussion and Conclusion 


With different scientific missions, some experimental spacecrafts have to 
operate with various kinds of slew motion for the purpose to perform its 
scientific experiments. In this study, the AXAF-S spacecraft has been given as 
an example to show a spacecraft operated with slew motion. The instability of 
the liquid-vapor interface surface of the fluid systems can be induced by the 
presence of gravity gradient and jitter accelerations associated with slew motion 
of the spacecraft. These instabilities originated from sloshing dynamics can 
cause various problems of spacecraft control systems. Sometimes, these sloshing 
dynamics problems can even deteriorate the quality of the normal operation of the 
spacecraft. It is vitally important to understand fully the characteristics of 
the spacecraft sloshing dynamics before one can assure the high quality operation 
of scientific spacecraft. 

In this report, the generalized mathematical formulation of sloshing 
dynamics for partially filled liquid in the dewar container driven by the both 
the gravity gradient and jitter accelerations associated with spinning and/or 
slew motions are investigated. Explicit mathematical expressions of both the 
gravity gradient and jitter acceleration associated with various types of 
spinning and/or slew motions are derived. The formulations of fluid moment and 
angular moment fluctuations in fluid profiles induced by sloshing dynamics, 
together with fluid stress and moment fluctuations exerted on the spacecraft 
dewar container have also been derived. 

Two examples of the numerical simulation of sloshing dynamics driven by 
orbital accelerations, applicable to the AXAF-S spacecraft associated with slew 
motion, have been carried out. These examples demonstrate (1) how great the 
degrees of sloshing dynamics for imbalance profiles have been developed; (2) how 
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great the amplitudes and fluctuations of angular momentum and moment have been 
induced; and (3) how great the forces and moment fluctuations exerted on the 
dewar wall of the container have been excited by the effects of gravity gradient 
and jitter accelerations associated with spinning and/or slew motions. The 
generalized mathematical formulations shown in this study can provide a useful 
tool for the development of effective control techniques to assure the high 
quality operation of the spacecraft to achieve the final goal of scientific 
missions . 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9 


Time -dependent variation of angular position of spacecraft slew motion 
which operates at 90° angular displacement in 10 minutes. 

Time -dependent variation of angular velocity of spacecraft slew motion 
which operates at 90° angular displacement in 10 minutes. 

Time -dependent variation of angular acceleration of spacecraft slew 
motion which operates at 90° angular displacement in 10 minutes. 
Three-dimensional cartesian coordinate (x" , y" , z") used for 

spacecraft slew motion with axes at mass center of spacecraft, and 
coordinates (x, y, z) used in fluid mechanics computation. 

AXAF-S spacecraft coordinate systems with azimuth angle from 

spacecraft mass center to the center of the Earth. Coordinate (x” , 
y" , z") for slew motion and coordinate (x , y, z) for fluid mechanics 
computation. 

Geometry of the GP-B dewar container with the coordinate system 
perpendicular and tangential to the container wall. (A) Geometry in 
r-z plane and (B) Geometry in r-0 plane. 

Geometry of the GP-B dewar container with baffle-boards and their 
locations. (A) Geometry in r-z plane and (B) Geometry in r-d plane. 
Distribution of grid points in the (A) Radial-axial plane, and (B) the 
radial -circumferential plane of the cylindrical coordinates for the 
AXAF-S dewar tank. 

Time variation of AXAF-S spacecraft gravity gradient acceleration 
acting on fluid mass located at (r, 0, z) - (12 cm, ir/2 , 3 cm) for 90° 
slew motion in 10 min. along the y"-axis and orbital period of 97.6 
min. (A) Along x-direction; (B) Along y-direction; (C) Along z- 
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direction. 


Figure 10 


Figure 11 


Figure 12 


Figure 13 


Figure 14 


Initial profiles of liquid-vapor interface for AXAF-S spacecraft of 
dewar tank under background gravity of 10~ 7 go and direction of 
background gravity at ^ — 0°. (A) In r-z plane at 6 - 0° and 180° , 
(B) In r-0 plane at $ - 90° and 270°, (C) In r -6 plane at z - 93.9 cm, 
and (D) Three-dimensional liquid-vapor interface profile. 

Time sequence evolution of AXAF-S spacecraft three-dimensional liquid- 
vapor interface oscillations for dewar driven by gravity gradient 
acceleration associated with slew motion in the y"-axis. 97.6 min 
orbital period and 90° slew motion in 10 min are applied to the 
spacecraft operation. 

Time sequence evolution of AXAF-S spacecraft liquid-vapor interface 
oscillations for dewar tank in r-z plane at 6 — 0° and 180°, driven 
by gravity gradient acceleration associated with slew motion in the 
y"-axis. 97.6 min orbital period and 90° slew motion in 10 min are 
applied to the spacecraft operation. 

Time sequence evolution of AXAF-S the spacecraft liquid-vapor 
interface oscillations for rotating dewar in r-z plane at 6 — 90° and 
270°, driven by gravity gradient acceleration associated with slew 
motion in the y"-axis. 97.6 min orbital period and 90° slew motion 
in 10 min are applied to the spacecraft operation. 

Time sequence evolution of the AXAF-S spacecraft liquid-vapor 
interface oscillations for rotating dewar in r-0 plane at height z - 
95.9 cm, driven by gravity gradient acceleration associated with slew 
motion in the y"-axis. 97.6 min orbital period and 90° slew motion 
in 10 min are applied to the spacecraft operation. 
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Figure 15 


Figure 16 


Figure 17 


Figure 18 


Figure 19 


Figure 20 


Time sequences of the AXAF-S spacecraft bubble mass center 
fluctuations due to the sloshing dynamics driven by gravity gradient 
acceleration associated with slew motion in the y-axis. 

Time variation of AXAF-S spacecraft jitter acceleration associated 
with slew motion acting on fluid mass located at (r, $ t z) - (12 cm, 
*/2, 3 cm) for 90° slew motion in 10 min. along the y"-axis and 

orbital period of 97.6 min. (A) Along x-direction; (B) Along y- 
direction; (C) Along z-direction. 

Time sequence evolution of AXAF-S spacecraft three-dimensional liquid- 
vapor interface oscillations for dewar driven by jitter acceleration 
associated with slew motion in the y" -axis. 90° slew motion in 10 min 
is applied to the spacecraft operation. 

Time sequence evolution of AXAF-S spacecraft liquid-vapor interface 
oscillations for dewar tank in r-z plane at $ - 0° and 180°, driven 
by jitter acceleration associated with slew motion in the y"-axis. 
90° slew motion in 10 min is applied to the spacecraft operation. 
Time sequence evolution of AXAF-S the spacecraft liquid-vapor 

interface oscillations for rotating dewar in r-z plane at 9 - 90° and 
270°, driven by gravity jitter acceleration associated with slew 
motion in the y"-axis. 90° slew motion in 10 min is applied to the 
spacecraft operation. 

Time sequence evolution of the AXAF-S spacecraft liquid-vapor 

interface oscillations for rotating dewar in r-0 plane at height z 
—95.9 cm, driven by gravity jitter acceleration associated with slew 
motion in the y"-axis. 90° slew motion in 10 min is applied to the 
spacecraft operation. 
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Figure 21 


Figure 22 


Figure 23 


Figure 24 


Figure 25 


Figure 26 


Figure 27 


Time sequences of the AXAF-S spacecraft bubble mass center 
fluctuations due to the sloshing dynamics driven by jitter 
acceleration associated with slew motion in the y-axis. 90° slew 
motion in 10 min is applied to the spacecraft operation. 

Time sequence of the AXAF-S spacecraft angular momentum fluctuations 
Hx, Hy, H z ) due to the sloshing dynamics driven by gravity gradient 
acceleration associated with slew motion in the y-axis. 

Time sequences of the AXAF-S spacecraft fluid moment fluctuations due 
to sloshing dynamics driven by gravity gradient acceleration 
associated with slew motion in the y-axis. (A) For moment along x- 
direction; (B) For moment along y-direction; and (C) For moment along 
z-direction. 

Time sequence of the AXAF-S spacecraft angular momentum fluctuations 
H*, Hy, H z ) due to the sloshing dynamics driven by gravity jitter 
acceleration associated with slew motion in the y-axis. 90° slew 
motion in 10 minutes is applied to spacecraft operation. 

Time sequences of the AXAF-S spacecraft fluid moment fluctuations due 
to sloshing dynamics driven by gravity jitter acceleration associated 
with slew motion in the y-axis. 90° slew motion in 10 minutes is 
applied to spacecraft motion. (A) For moment along x-direction; (B) 
For moment along y-direction; and (C) For moment along z-direction. 
Time sequence of the fluctuations of fluid stress forces exerted on 
the AXAF-S dewar container due to the sloshing dynamics driven by 
gravity gradient acceleration associated with slew motion. (A) Along 
x-direction; (B) Along y-direction; and (C) Along z-direction. 

Time sequence of the fluctuations of fluid stress moment exerted on 
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Figure 28 


Figure 29 


Figure 30 


Figure 31 


the AXAF-S dewar container due to the sloshing dynamics driven by 
gravity gradient acceleration associated with slew motion in the y- 
axis (A) Along x-direction; (B) Along y-direction; and (C) Along z- 
direction. 

Time sequence of the fluctuations of moment arm of viscous stress 
moment exerted on the AXAF-S dewar container due to the sloshing 
dynamics driven by gravity gradient acceleration associated with slew 
motion. (A) Along x-direction (pitching moment) ; (B) Along y-direction 
(yawing moment) ; and (C) Along z-direction (rolling moment) . 

Time sequence of the fluctuations of fluid stress forces exerted on 
the AXAF-S dewar container due to the sloshing dynamics driven by 
gravity jitter acceleration associated with slew motiony. (A) Along 
x-direction; (B) Along y-direction; and (C) Along z-direction. 

Time sequence of the fluctuations of fluid stress moment exerted on 
the AXAF-S dewar container due to the sloshing dynamics driven by 
gravity jitter acceleration associated with slew motion in the y-axis 
(A) Along x-direction; (B) Along y-direction; and (C) Along z- 
direction. 

Time sequence of the fluctuations of moment arm of viscous stress 
moment exerted on the AXAF-S dewar container due to the sloshing 
dynamics driven by gravity jitter acceleration associated with slew 
motion. (A) Along x-direction (pitching moment) ; (B) Along y-direction 
(yawing moment); and (C) Along z-direction (rolling moment). 
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